The NNPDF1.2 parton set: implications for the LHC 
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Recently a new set of Parton Distribution Functions (NNPDF1.2) has been produced 
^ — 1 and released by the NNPDF Collaboration. The inclusion of dimuon data in the analy- 

| sis allows a determination of the strange content of the proton with faithful uncertainty 

estimation together with a precision determination of electroweak parameters. In this 
contribution, we discuss some of the implications of the NNPDF1.2 set, and in par- 
ticular of its uncertainty determination of the strange PDFs, for LHC phenomenology. 
First of all, we study the impact on the electroweak boson production cross-section, 
O |' with special attention to the a(Z)/a(W) ratio. Then we revisit the top pair production 

cross-section, and perform a comparison of partonic fluxes between various PDF sets. 
Finally, we discuss the potential of using associated production of W with a charm 
quark at the Tevatron and the LHC to constrain the proton strangeness. 

; 1 The NNPDF1.2 parton set 

7-H ! 

• The determination of the strange and antistrange quark distributions of the nucleon is con- 

siderably interesting from the phenomenological point of view. However, till very recently, 
f*"*- ■ the bulk of data included in parton determinations, namely neutral-current deep-inelastic 

scattering, had minimal sensitivity to flavour separation, and no sensitivity at all to the 
separation of quarks and antiquark contributions. As a consequence, in standard parton fits 
the strange and antistrange quark distributions were not determined directly: rather, they 
were assumed to be equal and proportional to the total light antiquark sea distribution. 
Due to the availability of the new deep-inelastic neutrino and anti-neutrino charm pro- 
5^ ' duction data, which is directly sensitive to the strange and antistrange parton distributions, 

independent parametrizations of the strange and antistrange distributions have been in- 
cluded in most recent parton fits. 

However, the standard method for determining parton distributions, based on fitting 
the parameters of a fixed functional form, is known to be hard to handle when the ex- 
periments are relatively unconstraining. An alternative approach to parton determination 
which overcomes this difficulty has been developed by the NNPDF Collaboration in a series 
of papers The method is based on the use of neural networks for par- 

ton parametrization, and a Monte Carlo method supplemented by a suitable training and 
stopping algorithm for the construction of the parton fit. In this approach, parton distribu- 
tions are given as a Monte Carlo sample which represents their probability distributions as 
inferred from the data. 
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a See also Ref. [7] for a series of benchmark comparisons between the NNPDF approach and the standard 
approach. 
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Recently, in Ref. [5J, a new parton set, NNPDF1.2, was constructed^ The addition of 
dimuon data to the global inclusive deep-inelastic scattering dataset, on which the NNPDF1.0 [3] 
parton set was based, allows a determination of the strange and antistrange distributions 
with faithful uncertainty estimation. 

In Fig.[T]we compare the total strange PDF, s + (x, Qq), with Q\ = 2 GeV 2 , as determined 
from the NNPDF1.2 analysis with the results obtained from other parton sets: CTEQ6.6 [8], 
where an independent parametrization for s + is also used, and NNPDF1.0 and CTEQ6.5 [5], 
where the strange PDFs are fixed to be proportional to the non-strange sea distributions. 
We observe that in the case of the NNPDF analysis the main effect of releasing the strange 
parametrization is a substantial enlargement of the error on the s + distribution obtained 
in the NNPDF 1.2 analysis, which remains anyway compatible with the one of the NNPDF 
1.0 set. This is not the case for the CTEQ analysis where it can be observed that the 
s + distributions from CTEQ 6.6 and CTEQ 6.5 present substantial differences both in 
normalizations and in shapes, especially in the intermediate and small-re region. 

We refer to Ref. [6] for further details on the NNPDF1.2 parton set and for a detailed 
description of its implications in the precision determination of electroweak parameters, 
such as the CKM matrix elements \V CS \ and \V c d\ and the QCD corrections to the Paschos- 
Wolfcnstcin relation. 

Our purpose in this contribution is to briefly examine some of the phenomenological 
implications of the NNPDF1.2 parton set for LHC physics. In particular, we assess the 
impact on standard candles of the increased uncertainty in the strange PDF as compared 
to previous determinations. 

When comparing results computed using NNPDF sets with the ones obtained using 
other parton sets, like CTEQ 6.5/6.6, it is important to keep in mind that central values 
are affected by a systematic uncertainty due to an approximate treatment of heavy flavour 
contributions. On the other hand, the size of the PDF uncertainties on the same predictions 
should not be affected by the heavy flavour contribution. 

2 Impact of strangeness on LHC standard candles 

We begin by comparing results for the total cross-sections for weak bosons production at the 
LHC (cr (Vy ± ) and a (Z)), obtained using the NNPDF 1.2 set with the ones obtained using 
CTEQ6.6, NNPDF1.0, CTEQ6.5 and CTEQ6.1 [ID] sets. The computation of the cross 
sections has been performed including NLO QCD corrections using the MCFJV0 program. 
Results are collected in Tabled] 

It is interesting to notice that the size of the PDF uncertainties on these observables 
is very similar among all the different parton sets. This reinforces the observation already 
done in Ref. [8] that the uncertainty on total cross sections for weak boson production is 
rather insensitive to the uncertainty on the strange distribution. 

We turn then to a study of the sensitivity of the Z/W ratio rzw = a z/(^w+ + <J w-) 
to the uncertainty in the strange distribution. This observable is particularly interesting 
given that PDF uncertainties are greatly reduced when considering the ratio of two cross- 
sections, and thus provides an excellent candidate for a measurement of the LHC luminosity. 
However, in Ref. [8] it was shown that, although the impact of the uncertainties on the W 

b The NNPDF 1.2 set are available in LHAPDF starting from version 5.7.1. 

c Note that our results for CTEQ 6.6 do not coincide with the results presented in Ref. [8] because of the 
differences in the codes used for the computation of physical observables. 
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Figure 1: Comparison of the s + (x,Ql) PDF from various analyses. Upper left plot: NNPDF1.2 
vs. CTEQ6.6, where in both cases s is determined from experimental data. Upper right plot: 
NNPDF1.2, with fitted s + vs. NNPDFf.O, with s+ fixed by flavour assumptions. Lower left plot: 
CTEQ6.6, with fitted s + vs. CTEQ6.5, with s + fixed by flavour assumptions. Lower right plot: 
NNPDFf.O vs. CTEQ6.5, where in both cases s + is fixed by flavour assumptions. 

and Z cross sections due to the strange PDF arc rather small on their own, the ratio rzw has 
a greater sensitivity to the strange uncertainty in the region 0.01 < x < 0.05 because PDF 
uncertainties do not completely cancel in the ratio. Therefore, rzw is potentially affected 
by the larger uncertainty of strange PDF found in the NNPDF1.2 analysis with respect to 
other parton densities determinations. 

The NNPDF 1.2 prediction for the ratio r Z w and the correlation^ p [a(Z),a(W ± )j are 
given in Tab. [H together with the results obtained using the NNPDF 1.0, CTEQ6.5 and 
CTEQ 6.6 sets. In the Hessian approach, which is used for the CTEQ sets, the correlation 
between the two observables considered, o^W*) and cr(Z), is computed using the method 
described in [5]. In the Monte Carlo approach, used for computing PDF uncertainties for 
the NNPDF sets, the corresponding expression is given, as described in Rcf. g], by 

(0-(^)^(W /± ))rrn ~ ( CT (^))rrn (^(W ± ))r Pn 

p [a{Z),a{W ± )} = rcp rcp rcp , (1) 

V>(^) 2 } rcp - (°{z))l^{o{w±y)^ - M^±)) r 2 cp 

where the averages are performed over the N Ycp replicas of the NNPDF sets. 

d Ref. [8] instead uses the notation cos ip to denote the correlation between two PDFs/obscrvablcs. 
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er(W+)Br (W + -> l + Ui) [nb] 


a(W-)Br [W~ -» i+i/j) [nb] 


a{Z")Br (Z" — [nb] 


NNPDF 1.0 


11.83 ± 0.26 


8.41 ± 0.20 


1.95 ± 0.04 


NNPDF 1.2 


11.99 ± 0.34 


8.47 ± 0.21 


1.97 ± 0.04 


CTEQ6.1 


11.65 ± 0.34 


8.56 ± 0.26 


1.93 ± 0.06 


CTEQ6.5 


12.54 ± 0.29 


9.19 ± 0.22 


2.08 ± 0.04 


CTEQ6.6 


12.41 ± 0.28 


9.11 ± 0.22 


2.07 ± 0.05 



Tabic 1: The LHC benchmark cross sections for NNPDF1.0 and CTEQ6.5 (with strangeness 
proportional to the non-strange sea) and NNPDF1.2 and CTEQ6.6 (with strangeness determined 
from the global analysis). The CTEQ6.1, with a ZM scheme for heavy quarks, is also shown 
for comparison. All numbers shown are for y/s=lA TeV. PDF uncertainties correspond to 68% 
confidence levels. 





NNPDF1.2 


NNPDF1.0 


CTEQ6.6 


CTEQ6.5 


rzw 


0.0961 ± 0.0005 


0.0965 ± 0.0003 


0.0964 ± 0.0004 


0.0957 ± 0.0002 


p\a(Z),a(W ± )\ 


0.976 


0.994 


0.983 


0.994 



Table 2: Comparison of the values for the ratio of the Z and W cross sections at the LHC as well 
as their correlation p [a(Z), a(W ± )] , Eq. [TJ computed with different PDF sets. Again, all results 
are computed for a center-of-mass energy of ^/s=14 TeV. 

In Fig. [2] we compare the oz-gw one sigma correlation ellipses for the NNPDF 1.2, 
NNPDF 1.0, CTEQ 6.6 and CTEQ 6.5 sets. We note that, despite the fact that the error 
band on the strange parton densities is in general much larger for the NNPDF 1.2 set than 
for CTEQ 6.6, the uncertainty on the ratio rzw is of the same size. This is a consequence 
of the fact that, as previously mentioned, this ratio is mostly correlated to the strange 
PDFs in a limited region of x, namely 0.01 < x < 0.05, where the NNPDF 1.2 and CTEQ6.6 
uncertainties on s + are roughly of the same size. This can be understood as a consequence of 
the fact that the NuTeV dimuon data, which constrains the strangeness in the two analyses, 
cover precisely this kinematical range. 

The situation could be different if we look at the the differential rapidity distribution. 
In order to check this, in Fig. [3] we show the rapidity distribution of the ratio rzw denned 
as 

dr Z w , , _ dcr z (y)/dy 
dy (y) ~ d^{y)/dy' W 

together with the associated PDF uncertainties. We observe a sizable increase in the PDF 
uncertainty for the NNPDF 1.2 result, when compared to results obtained with other sets, 
at forward rapidities. This is due to the increase of the NNPDF1.2 strange uncertainties at 
small-x, shown in Fig. [T] However in the central rapidity region, which provides the bulk 
of the contribution to rzw, the uncertainties of CTEQ6.6 and NNPDF1.2 turn out to be 
comparable, confirming the agreement of PDF uncertainties on the integrated ratio shown 
in Tabic H 

3 The ti total cross section and partonic fluxes 

The total ti cross section, a(ti), has been the subject of several recent studies [8l [TTl [L2] . 
partially motivated by the possibility of using this process as a standard candle at the LHC 
to measure the luminosity. To compare the predictions of the NNPDF1.2 set to ones of the 
other parton determinations, the total ti cross section is shown in Fig. |4] for various PDF 
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Figure 2: Comparison of the W and Z one sigma correlation ellipses at the LHC obtained from 
different fits: NNPDF 1.2 (green), NNPDF 1.0 (blue), CTEQ 6.6 (red) and CTEQ 6.5 (purple). 



sets. 

The computation has been performed with MCFM at NLO. No soft gluon resummation 
corrections, as done in Ref. [TTj . are included. All scales are set equal to m t = 172.5 GeV. 
Together with NNPDF1.2, we show the predictions from NNPDF 1.0 and various other sets 
from the CTEQ and MRST/MSTW collaborations. It is clear from Fig. gj comparing for 
example NNPDF 1.0 and 1.2, or CTEQ 6.5 and 6.6, that strangeness plays a rather minor 
role for this observable. 

In order to understand the results of Fig. [4j one should take into account the fact that 
the values of a s (M§ ) used each global PDF are different. This is especially important for 
hadronic cross sections in which the leading order process is a strong interaction process (like 
in ti or gg — ► H), since in these cases the cross section at leading order will be proportional 
to a 2 s . Moreover, in general there is a correlation between a s and the PDF uncertainties, 
mostly the gluon, as discussed for example in Refs. [T31Q3]. Note that, while CTEQ 6.6 
and NNPDF 1.2 use the value of as(M^) roughly along the lines of the most up-to-date 
global average (a s = 0.118 and a s = 0.119 respectively), MSTW08 Q3] determines it from 
the PDF analysis ending up with a s = 0.1202^°;°°^ at NLO. 

The use of different values of as has to be taken into account when comparing this 
kind of standard candle processes. For example, the ratio [a™ stw08 /a^ tGq66 ] , which appears 
in the Born top pair production cross-section, implies that, even for identical PDFs, the 
MSTW08 cross-section will be ~ 4% higher than the CTEQ6.6 one. Therefore, the use of 
different values of as could explain part of the discrepancy between various PDF sets. This 
emphasizes the importance of using a unique value of the strong coupling for comparison 
between benchmark observables, especially those proportional to as already at leading order. 

In order to compare the contributions from each PDF flavour to the total cross-section 
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The Z/W ratio 



The Z/W ratio 
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Figure 3: Left: the differential rapidity distribution of the ZW ratio, defined in Eq. [2] The very 
small PDF uncertainties are not shown. Right: the relative PDF uncertainties in rzw, as a function 
of the rapidity y. 



for tt production we follow the approach of Rcf. [15] . where the hadronic cross section is 
written as 



a (S, m? 



E 



— 



9 



4m| 



(3) 



that is, as a convolution between a partonic cross section cr^ and parton fluxes $y , defined 
as 



o Jo 



®ij (t,^ 2 ) = / dx\ \ dx 2 qi (xi,fi 2 ) q 3 (x 2 ,/i 2 )J (xix 2 - r) 



(4) 



which in turn are convolutions of parton distributions. Note that at the LHC, for \/~S = 14 
TeV, for a top quark mass of m t = 172.5 GeV one obtains p* ~ 6 • 10 -4 . 

The dominant contribution to top pair production at the LHC comes from the gluon- 
gluon and gluon-quark channels. In Fig.[5]we show the absolute fluxes for these two channels, 



as defined in Eq.|4]for /i 2 



for the three most recent PDF sets of the CTEQ, MSTW [16] 



and NNPDF Collaborations, in the kinematical region relevant for top pair production. The 
three PDF sets considered agree reasonably well. In order to obtain a more detailed picture 
of the comparison, in Fig. [5] we show the relative differences between the PDF fluxes for 
the two channels considered before, using the central CTEQ6.6 result as a reference. We 
observe a good agreement for the fluxes in the GG channel, as well as for the QG channel 
in most of the r range, with the exception of the middle region where a 2-rj discrepancy is 
found. The origin of such discrepancy might be related to differences in the treatment of 
heavy quark mass effects in NNPDF1.2 compared to the other sets. 

Finally, let us compare in more detail the size of PDF uncertainties in the parton fluxes. 
Following Ref. [T7], we compare the relative PDF uncertainties from the fluxes computed 
with NNPDF 1.2, CTEQ6.6 and MSTW08 as a function of S = xix 2 S. That is, we compute 
PDF uncertainties in $y(r = S/S,/j, 2 ) as a function of S. These parton fluxes have been 
computed at /j 2 = 10 4 GeV 2 , the typical scale for processes like W, Z or H production at 
the LHC, again assuming V~S = 14 TeV. 
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Top pair production at LHC, MCFM, m, = n R/F = 172.5 GeV 
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Figure 4: Comparison of the total ti cross section computed with recent PDF sets. The computation 
has been performed with the MCFM program at NLO for mt = /ljk = Hf = 172.5 GeV. All PDF 
uncertainties correspond to 68% confidence levels. 



We show the results of the comparison in Fig. [7] Note that only PDF uncertainties are 
shown, all central values are set to 1, unlike the case of Fig.[5j We observe a reasonable agree- 
ment in the size of PDF uncertainties in the intermediate S region, and sizable differences 
at smaller S. fn the intermediate S region, relevant for the production of massive objects, 
the PDF uncertainty in the QG channel turns out to be somewhat larger in NNPDF1.2 as 
compared to both CTEQ6.6 and MSTW08. This could be due to missing hadronic data 
in the former or to a parametrization bias in the latter: the upcoming NNPDF2.0 global 
analysis should settle this issue. 

4 Wc production and constraining strangeness at the LHC 

In the last section of this contribution, we focus on a LHC process which in principle could 
be used to measure the strange PDF which is presently poorly constrained in the region 
x < 10~ 2 [6], Wc associated production. The associated production of a vector boson and 
a charm quark at hadronic colliders is directly sensitive to the strange sea PDF and for this 
reason, in the past it has been proposed as a candidate for constraining the strange PDF at 
the TeVatron and the LHC colliders [15] , 

We revisit this proposal by comparing the NNPDF1.2 predictions for the total Wc cross 
section with recent measurements at the Tevatron |19j and giving predictions for the LHC. 

The dominant production channel for Wc production is gq — > Wc, with q a down-type 
quark. As in the case of neutrino dimuon production, the down and bottom initiated contri- 
butions are suppressed with respect to the strange one by the smallness of the corresponding 
CKM matrix elements, and therefore at LO, neglecting CKM mixing, the cross section is 
proportional to the strange PDF. 

Associated W and charm production looks therefore like a promising channel for pro- 
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Figure 5: The partonic fluxes, defined in Eq.|4j in the QG (left) and the GG (right) channel, for the 
NNPDF1.2, CTEQ6.6 and MSTW08 PDF sets. The rapid decrease at large r reflects the smallness 
of the PDFs at large x. Note that the partonic fluxes are evaluated at = mf . 




Figure 6: The relative differences between the PDF fluxes in the QG and GG channels between 
the various sets considered, in the kinematical region relevant for top pair production at the LHC 
As in Fig. [5j partonic fluxes are evaluated at fi 2 = mf. Only the kinematically allowed region 
t > T m i n = 4m 2 /S is shown. 



viding a direct constraint on the strange PDF at the energy scale of Mw, two orders of 
magnitude above the typical energy of the NuTeV dimuon data. 

This process has already been studied at the Tevatron. In particular, the CDF experi- 
ment has published a measurement of the Wc production cross section, obtained using ~ 1.8 
fb _1 of pp collisions at y/s = 1.96 TeV. A NLO prediction for this observable can be obtained 
easily using the MCFM gD] code. 

The result obtained by CDF for the measurement of the total cross section is 

owc{ptc > 20GeV, \r) c \ < 1.5) x BR(W -> lu) = 9.8 ± 3.2pb, (5) 

which is in agreement with the NLO prediction obtained using MCFM and the NNPDF 1.2 
set: 

<JWc(ptc > 20GeV, |?7 C | < 1.5) x BR(VF -► lu) = 10.11 ± 1.24(PDF)tg;^(scale) pb, (6) 
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Figure 7: Comparison between the size of PDF uncertainties in the partonic fluxes in the 
NNPDF1.2, CTEQ6.6 and MSTW08 PDF sets in the region 1CT 5 < x < 1 (above) and a zoom in 
the middle region (below), where PDF uncertainties are smallest. Note that in this case partonic 
fluxes are evaluated at fi 2 = 10 4 GeV 2 , as discussed in the text. 



where the first error is the one coming from PDF uncertainties and the second one is due to 
the variation of the renormalization and factorization scales in the pcrturbative computation. 
The expected precision of the experimental result, when extrapolated to the full Run II 
dataset (^6 — 7 fb _1 ), is ~ 15%, comparable to the present uncertainty on the theoretical 
prediction. The theoretical uncertainties at the Tcvatron is dominated by PDF uncertainties, 
with a sizable contribution from scale dependence. 

In order to investigate the possibility of using this very same channel as a strangeness 
constraint at the LHC, we also computed the Wc cross section for the LHC assuming a 
center of mass energy of 14 TeV and standard rapidity and transverse momentum cuts both 
for the charm quark and the leptons coming from the W decay. The result that we obtain 
is 

(?Wc(ptc > 20GeV/c, \r) c \ < 4.) x BR(1F Iv) = 631 ± 46 (PDF) ±H (scale) pb . (7) 

Our result shows that the uncertainty on the theoretical prediction due to scale variations, 
and thus to higher order corrections, is comparable to the uncertainty due to the strange 
PDFs. The result seems to suggest that it might be difficult to use this process to constrain 
PDFs. It might instead be useful to look at differential distributions, like the W or c- 
jet rapidity or transverse momentum distributions. In any case, unless scale uncertainties 
can be reduced by higher order computations, the theoretical error limits the usefulness of 
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Wc production as a constraint for the strangeness distributions, regardless the accuracy of 
present and future experimental measurements of this channel. 

5 Outlook 

In this contribution we reviewed some phcnomenological implications of the NNPDF1.2 
parton set for LHC physics, with emphasis on those observables more directly sensitive 
to the strange PDFs. A more detailed study of the phcnomenological implications of the 
NNPDF approach to LHC observables is however postponed until the release of the upcoming 
NNPDF2.0 global parton analysis, which includes all relevant hadronic data like Drcll-Yan 
pair production, vector boson production and inclusive jet production. 
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